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Abstract 



General quantum mechanical states of neutrinos produced by mechanisms outside the Standard 
Model are discussed. The neutrino state is described by the Maki-Nakagawa-Sakata-Pontecorvo 
unitary mixing matrix only in the case of relativistic neutrinos and Standard Model left-handed 
charge current interaction. The problem of Wigner spin rotation caused by Lorentz transformation 
from the rest production frame to the laboratory frame is considered. Moreover, the mixture of the 
neutrino states as a function of their energy and parameters from the extension of the Standard 
Model are investigated. Two sources of mixture, the appearance of subdominant helicity states, and 
mass mixing with several different mixing matrices are studied. 
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1 Introduction 



The neutrino oscillation phenomenon is well established. In several different experiments, the transition 
between different neutrino flavours is observed There is agreement that oscillation is possible when 
different mass neutrino states add coherently [2] in their production process. Such a combination of eigen- 
mass neutrino states was defined many years ago by Maki, Nakagawa and Sakata using the previous 
idea of neutrino mixing invented by Pontecorvo|4i. Now, there is agreement concerning the definition of 
so-called neutrino fiavour states for relativistic neutrinos. The state of three massive relativistic neutrinos 
produced in the process involving a charged flavour a lepton is given by [S] 



where U is the unitary mixing matrix known as the Maki-Nakagawa-Sakata-Pontecorvo matrix. Such 
a state, which describes neutrinos with negative helicity, is understandable. Neutrinos are produced in 
charged lepton processes, which are described by the chiral left - handed interaction; the only one existing 
in the Standard Model (SM). Coupling of the charged la lepton with a massive i/i neutrino and W boson 
is given by: 



Such an interaction produces a relativistic neutrino with negative helicity (antineutrino with positive 
helicity), and its propagation in a vacuum or in matter does not reverse spin direction (neutral current 
interaction is also left-handed). In such circumstances, it is natural to describe neutrinos by a pure 
quantum mechanical state Eq.([T]). So far, only relativistic neutrinos are detected (the detection energy 
threshold is practically larger then 100 keV [6]), and the pure production state Eq.(IT]) works well. However, 
there are several reasons to pose the following question about general neutrino states: What is the neutrino 
state if it is not necessarily relativistic, and/or its interaction is more general than that given by Eq.(l2])? 
The answer for such a question is important as attempts are ongoing to define different neutrino states 
in a non-relativistic region (see e.g. [7J ). Besides, there is a chance that future more precise neutrino 
oscillation experiments will open a window for New Physics (NP) and the knowledge of the initial and final 
neutrino states could be important. The attempts to define general neutrino states have been discussed 
many times in the literature (see e.g. [8^). The neutrino production and detection states that depend on 
the specific physical process were considered in e.g. 9 . To our knowledge, the most general oscillation 
and detection neutrino states were studied by Giunti lOJ (see also [Hi ). However, even in this work, 
only pure states in the relativistic case, without the spin effect of particles which produces neutrinos, are 
finally discussed. 

The aim of our study is to construct a general neutrino density matrix using the helicity - mass basis 
using the amplitude for the production process in the CM system or rest system of a decaying object 
(e.g. muon in the neutrino factory or nucleus in the beta beam experiments) using a full wave packet 
description of particles. In order to find how physics beyond the SM modify neutrino production states, 
we consider the general neutrino interaction. Next, we transform the density matrix to the frame of the 
detector - the laboratory frame. Such a general Lorentz transformation causes Wigner rotation for spin 
states [12]. We show that, even if such a Wigner rotation can be large, its effect on neutrino helicity states 
is small and without practical meaning. Then we have the general density matrix describing a neutrino 
beam, which oscillates as it passes between the production site and the detector. In this way, we will be 
able, although it is not subject of this paper, to consider more precisely any possible infiuence of NP not 
only on oscillation phenomena as usual ([13]), but also on the production and, in a similar way, detection 
of neutrino. In such a case, all quantum mechanical properties of neutrino oscillation phenomena (such 
as the problem of factorization, oscillation and coherence length) could be considered in the more precise 
way. 

We ask questions about when the neutrino state is pure and when it must be consider as mixed, how 
this property depends on neutrino energy, and the mechanism of the production process. We show that, 
for non-relativistic neutrinos, their state is always mixed. This is a very natural result, for neutrinos with 
an energy comparable to their masses, even a left-handed interaction (Eq.([2])) produces neutrinos with 
both helicities. Such a case must be described by a mixed state. But, as we will see, the appearance 
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of the subdominant neutrino helicity states is not the only reason for the mixture of states. If the NP 
introduces different mixing matrices for different kinds of interactions (for vector and scalar couplings, 
as examples), quantum mixing of states also appears . Because of that, any attempt to treat a flavour 
neutrino as one object in Fock space (see e.g. [7]), even if mathematically consistent, creates a problem 
with the description of the real production process. Neutrino states are always mixed if the production 
mechanism is more complicated in comparison with the SM left-handed interaction Eq.([2]). The neutrino 
density matrix also depends on the other properties of the production process; for example, the scattering 
angle, polarizations of other accompanying particles, and their momentum uncertainties as well as on the 
neutrinos themselves, the mass of the lightest neutrino, or their mass hierarchy. 

In this paper, we discuss only the first stage of the oscillation process - the preparation of neutrino 
states - leaving the problem of oscillation, together with the interaction inside the medium and the 
problem of their detection, for a future paper. In the next Chapter, we describe the neutrino density 
matrix for the selected production process. As the base process, we take the neutrino production where 
charged flavour leptons a scatter on nuclei A (la + A ^ Ui + B). Amplitudes are calculated for the 
general interaction and the normalized {Tr{p) — 1) density matrix in the CM system is found. Then, 
the problem of the Lorentz transformation and Wigner rotation is discussed. The full analytical and 
numerical properties of the neutrino density matrix are analyzed. To distinguish between mixed and 
pure states, we calculate Tr{p'^) as a function of neutrino energy and the parameters which describe the 
SM extensions. We also discuss how fast the neutrino states become pure as the energy increases. We 
study density matrices for two subsystems: the first which describes both neutrinos helicity independent 
of their mass, and the second describing the mixture of three neutrinos, independent of their spin state. 
Finally, Chapter 3 gives our final comments and future prospectives. 



2 Neutrino states in the production processes 

Let us assume that the beam of neutrinos with mass i, Vi is produced by charged leptons (a — e, /i) which 
scatter on nuclei A. 

la + A^v, + B. (3) 

We do not measure separate neutrinos and their masses; in fact, we do not directly measure neutrinos at 
all. Quantum neutrino states can be deduced from the measured states of all other particles accompanying 
the neutrino. So, having information about particles A, B and, especially, about charged lepton in 
the process Eq.([3]), gives us the ability to define neutrino states. To define the spatial or momentum 
distribution for neutrinos, we have to assume that all particles in the discussed process do not have a 
precisely determined momentum, so they must be defined by wave packets and not by plane waves (as 
is usually the case in other scattering processes). In order to discuss coherence effects in the detection 
process, we have to define the momentum distribution of the mass-helicity density matrix. Such a density 
matrix can be constructed from the amplitude for the process Q whose interaction mechanism is under 
investigation. 

We assume that our interaction is described by an effective model Lagrangian where both the vector 
and scalar couplings with the left and right chirality interaction are present. This sort of interaction 
appears, for example, in the Left - Right symmetric model (see e.g. [2]). We assume that: 



ccc= - ^ J". ^ \y.^.w^^- "f^VL^^: + 7^(1 + 75)4£/„t] u w+ 

2\/2sin0w i'~ 
+ 5] [ (1 - j,)ijlV^* + (1 + J5)vrV^*] L H+ 

+ Y,^[i^-l5)rLW^d* + {^ + J5)rRW^;]dH+]+h.c, (4) 



where eL,R, Vl,r and tl,r are parameters that are taken to define the scale and differ slightly from their 
lySM values. 
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In the neutrino sector, we introduce the mixing matrices separately for chiral-left and chiral-right 
parts as well as for vector and scalar neutrino interactions. As a consequence, phenomena such as flavour 
neutrino states, which were explicitly defined, do not exist any more, even for relativistic neutrinos. 
Mixing in the left-handed vector interaction (?7^) and left-handed scalar one {V^) can be different and 
an unambiguous definition of flavour neutrino state is impossible. We see that the determination of the 
flavour neutrino state given by the Pontecorvo-Maki-Nakagawa-Sakata matrix is purely accidental and 
valid only in the SM and relativistic neutrino interaction. In such circumstances, we cannot consider 
flavour neutrinos as elementary quanta of a field theory. 

In the lowest order, two Feynman diagrams describe the neutrino production process Eq.Q 




Figure 1: The charged boson(s) and charged Higgs particle(s) H'^ exchange diagrams, which describes the 
neutrino production process. 



To define the neutrino state we have to calculate the amplitude for the process Eq. ([3]). Three 
measured particles la, A and B are in the states given by the wave packets (a = /, A, B) 

IV'a,Aa)= d^Pa'^aiPa,PaO,Crap)\Pa,K): (5) 



where \pa, ^a) is the particle state with definite momentum pa and helicity Xa, and 5'a(Pa,PaO; crop) are 
the momentum distribution wave functions with a central momentum pa and momentum uncertainty aap- 
For A narrow momentum distribution {aap <C \Pa\), the physical results are independent of the specific 
shape of these functions, and we take them in the Gaussian form; 

'i!a{Pa,PaO,Crap) ^ — —3-6 """p ■ (6) 

(27rcr2p)4 

With such wave packet states, we calculate the amplitude for the process ([3]). After integration over the 
particles momenta we obtain: 

{f\S\^) - - V,^ I , / d^xe-^(^-^-^'e^(P-P->^e-^ ""C^"' Af(p,A,...), (7) 

where the average values of neutrino energy and momentum are determined from the appropriate con- 
servation rules (we use the same notation as in Ref ([IH])): 

Er = E{pio) + E{pAo) - Eip^„), 

and 

Pfl ^Pio+PAo -Pbq- 

For the precise energy and momentum determination we have to remember that, for a given total en- 
ergy in the process ([3]), the final particle momenta and energies depend on the neutrino masses. The 
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neutrino spatial uncertainty af/x and average velocity vr are obtained from the accompanying particles 

1 

2o-„ 



uncertainties aax — 2^— and their velocities Va — {^)p=pa 



And, finally, 
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The amplitudes M(p, A, ...) = Af (p. A; A;, A^, As) are the full helicity amplitudes for the production 
process Eq.Q, which is calculated in the CM frame for average values of all particle momenta except 
for the neutrino. Using the saddle-point approximation, the integrations over space and time can be 
done, giving the final amplitude Af{p, A, Xb', A;, Aa) for the production of the Vi neutrino with mass m^, 
momentum p, and helicity A in the process initiated by flavour a lepton: 

Af{p, A, As; A,, A^) = E) M^{p, A; A;, A^, As), (8) 

where ^ ^ 



' Rp 

and 



fUp,E) = e ■ 



7^ I2 ,-r2 \ 



where we have openly indicated that both factors depend on the neutrino mass. Once more, momentum 
uncertainty anp = appears in the above formulas. The second function fai describes the momentum 
distribution for the produced neutrinos. 

The amplitude Eq.® describes all properties of the produced neutrinos. Their quantum mechanical 
state is determined by the density matrix defined as: 

z. A; p',/,a') = -^ ^ A':{p,\Xb-AiAa)AT{p' ,\,\b;\iAa), (9) 

where we assume that, in the production process ([3]), the initial particles are not polarized and the final 
polarization of the B nuclei is not measured. For the normalization factor iVo,, given by 

3 

iVa=EE E d^P\A^{P,\^B-AiAA)\\ (10) 

A=±l i=l \i,\a,\b 

the density matrix is properly normalized: 

3 

TriQ'^) ^ h^'P e"(p,A,z;p,A,z) = 1. (11) 

A=±l i=l 

Eq. ([HI) describes the neutrino state in the CM frame. Such a frame is equivalent to the rest frame of a 
decaying muon in the neutrino factory case or the decaying nuclei in the beta beam experiments. Usually, 
neutrinos are not produced in such a frame. They are produced in the laboratory frame (LAB), where 
decaying muons or nuclei are moving very fast. The density matrix in such a frame is interesting for us. 
To reach the LAB frame, we have to make a Lorentz transformation. As helicity is not a full Lorcntz 
scalar, we can expect that the density matrix in the LAB will be rotated (helicity Wigner rotation). Each 
neutrino state \p, A, i) with mass m^, momentum p and helicity A after a general Lorentz transformation 
A, does not necessary point along the neutrino momentum, but transforms as: 

U[AUs,X)=Y,^ax(R'^{^,'5,-^)RwR{v,d,-v))\Ap,s,a), (12) 
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where angles ip and 9 {ip and 6) are the spherical angles of momenta before (p) and after (Ap) which 
carry the Lorentz transformation. R\y is the normal Wigner rotation: 

Ry^, = L-^{Kp)KL{p). (13) 

Hence, each element of the density matrix transforms as: 

A; p ,z ,A ) -> Q"{Kp, i, 77; Ap ,i ,t]) = (14) 

^Dy,^(i^;VJ^'"(P, h A; p',t',X)'Dl('*{R'^.), 
aa' 

where is the helicity Wigner rotation: 

= R^i^,o, -ip)RwR{^, 0, -(^), 

and depends on neutrino mass as the neutrino momentum depends on it. The matrix 'D^J'^{R^^^) can be 
found. If we parameterize: 

a 

'b' a 



D^/'«.) = ( t t ] , (15) 



then 
where 



\b\=- 



b=\b\e-^^, 
S{p + CE cos(6l)) - CE sm{6) 



' 2[{E +p cos{0)){E + Cp + p cos{9) +CE cos{e) + SE sin{e))] 
and 

p = p/37, E = mi + i?7, 

with 

- ^ 7(i^/? + pcos(.)) ^.yiTT^. 
^"f^{EP+pcos{9))^ +p2 sin2(0) 

Here m^, and p are the neutrino mass, energy and momentum, 7 and /3 are parameters of the Lorentz 
transformation in the "z" direction, and 0, ip are spherical angles of the neutrino momentum. It is worth 
noticing that, for normal (not helicity) Wigner rotation Rw, G — > cosO and S — > sinO. The off-diagonal 
element of the D^/^ matrix Eq. (fT51) b decides how big the effect of the Wigner rotation is. In Fig. 
2, on the left side, we have plotted the |6| element as a function of neutrino 7 factor (7^ ~ ■^) for 
various scattering angles 9cm for the Lorentz transformation along the "z" axis with 7 = 500. We see 
that the effect is large but only for very small neutrino energies {E Ki rrii) and large scattering angles. 
Unfortunately, as we see from right side of the picture, such neutrinos also have small energies in LAB 
frame and are not measured. The effect of the helicity flip for experimentally measurable neutrinos with 
a enough large energy [Elab ^ 100 MeV) is completely negligible. We see that, for practical purposes, 
we can safely take it that Lorentz transformation does not change the helicity structure of the density 
matrix and 

f?^^^(Ap, i, A; Ap',z',A')= e^*^(p, z. A; p',z',A'). (16) 

We can now check the quantum structure of the neutrino state described by Eq.Q. First, we would like 
to answer how good the MNSP assumption (Eq.(IT])) is. In order to do that, we calculate Tr[{Q°'y'\: 



Tr[{Q''f\ ^ J2 J'^'P dh'g'^iP. h A; p ',^',A')^?"(p ',A'; p, 



i. A), (17) 



Ai A 

and check how it differs from unity. 
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= 1 eV, Y - 500 




Figure 2: On the left - the off-diagonal element |6| of the helicity Wigner rotation matrix as a function of the 
neutrino's 7^ — factor for five different neutrino scattering angle ^ = f , f , ^ and ^ in CM. The Lorentz 
transformation is described by the 7 — 500 factor. On the right side, the neutrino momentum in LAB system for 
the same scattering angle in the CM frame is shown. 



To be more precise and to not complicate our consideration with the structures of particles A and B, 
which are not important for our purpose, we take the u and d quarks with masses m„ = 2.25 MeV and 
nid — 5.0 MeV as particles A and B, respectively. For lepton a we take an electron with rue = 0.511 MeV. 
All Particles momenta uncertainties are equal and Cep = crup = crdp = 10 eV. The results do not depend 
crucially on the values of these uncertainties. In Fig. 3, we plot 1 — Tr'[(g")^] as a function of the kinetic 
neutrino energy for two values of the lightest neutrino mass, mi — eV and toi = 1 eV in two cases, 
when only Left-Handed chiral interactions are present (only 7^ 0), and when Left- and Right-current 
exist (el and ^ 0). The masses of the heavier neutrinos and the vector Left-Handed mixing (Ul) 
are taken from the oscillation data [TS]. We assume a normal mass hierarchy with ttiq = 8.0 x 10^^ eV^, 
m^(„ = 2.5 X 10~^eT^^ and mixing angles, sm^0i2 — 0.31, sm^023 = 0.44, sin^Oi^ = 0.009 and 5cp = 
(central values are taken). The mixing in the Right-Handed vector {Un) and scalar interactions {Vl^r) we 
take to be the same as in the L-H case [Vl = Vr = Ur = Ul)- As we will see, by taking the same mixing 
matrices, we ignore the QM mixing caused by different mass states, and assume possible mixing effects 
are caused only by different helicity states. If the production process is described by the SM interaction 
(only L-H vector interaction), the state is strongly mixed for very small neutrino energies (continuous 
line and dotted line). The effect of the mixture depends on the lightest neutrino mass (as we see from 
the Fig. 3) and depends also slightly on the mass hierarchy and chosen particle momentum uncertainties 
(not shown). 

If the neutrinos energies increase, their state becomes pure. For relativistic neutrinos where masses 
and momentum uncertainties can be neglected, the density matrix takes the form: 

^,"(+1, z; +1, fc) = ^ (M 4 U^^U^, + A^. r,l V^*V^, + ertjl {U^*V^k + Vt*U^,)) , (18) 

f?"(+l,i;-l,A;) = £)"(-l,fc;+l,z)* = (19) 

a. 

e"(-l, z; -1, fc) = ^ [M el U^Utk + A^.^ rjj, V^V^ + A,,,, slVr {U^tV^^ + V^CUD) , (20) 

and 

= A,i el + iql + A,2^ el + ri + (21) 
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Figure 3: 1 — Tr[p^] as a function of the neutrino kinetic energy for the two lightest neutrino mass mi — eV 
and mi = 1 eV with only Left [el = l,eR = 0.) and the Left- and Right-currents (el = l,eR = 0.01). 



where the A factors (not given here) depend on the neutrino energy and scattering angles. We see that, 
in the SM where only 7^ 0: 

e"(-l,z;-l,fc) = C/i^*C/ffe, (22) 

and all other elements of density matrix are equal zero. Such a density matrix describes pure states 
equivalent to the MNSP mixture Eq. ([T]). Relativistic neutrinos in the SM model are produced 
in the pure states described by the MNSP mixture. The situation changes if any NP participates 
in the neutrino production process. Then, the states are not pure - even for relativistic neutrinos. A 
departure from the pure state depends on two properties of the system: the helicity mixing and the mass 
mixing. The first is obvious - neutrinos in two helicity states are produced which, in a natural way, gives 
the mixed state. The second effect is more mysterious and appears only if the mixing matrices are not 
the same. In order to investigate both effects, we can construct separately two density matrices for two 
subsystems - the first for the helicity subsystem and the second for the mass. For the helicity subsystem, 
the density matrix is given by: 

eA:r" = E^^"(^'^;^'*)' (23) 

1=1 

and, for the mass subsystem: 

<r^^= E e"(A,»;A,fc). (24) 

A=±l 

With reference to the helicity subsystem, 1 — rr[(p"^''P™)^] is different then zero even if all mixing matrices 
are equal, and it is a second order effect in the NP parameters: 

l-rr[(e"'^P")2]«0(4,4,ryi). 

In Fig. 3, we see (dashed line (toi = QeV) and dashed - doted line (mi = leV) that above the non- 
relativistic region, where even the L-H vector current produces neutrinos with h=+l/2, the state mixture 
is almost constant and of the order of £|j « 0.0001. Generally, the mixture of the neutrino states is, in 
a good approximation, proportional to the probability squared that neutrinos with negative helicity are 
produced. 
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Tr[(e")2] « (rr[g"P_])2, (25) 
where the P_ is the projection operator on the negative heUcity state, fc = ^A.-if^r/.-if^i.fe- In 




Figure 4: l-rr[(£i"'""'^'')^] as a function of the kinetic neutrino energy for various mixing matrices in the vector 

and scalar Lagrangian. For the parameter which defines New Physics, we take er — 0.01, 771, = rjn = 0.02. 
All matrices are parameterized in the same way as the standard Ul matrix. a) sin^9\2 = 0.87, sin^O^s = 
0.09, sin^9i3 = 0.23, Scp = 0.5 for Ur and Vl = Vr = Ul; b) Vl = Vr as Ur in the case a) and Ur = Ul; c) 
Vl = Vr the same as in b) and sm^Si2 — 0, siv?62-A ~ 0.09, siri^Ois — 0.71, Scp = 0.8 for Ur; d) Vl = Vr as in 
the case b) and sm^6'i2 = 0.01, sin^Oiz = 0.09, sin^dw = 0.71, 5cp = 0.8 for Ur. 

Fig. 4, we depicted the effect of neutrino state mixtures caused by the mass mixing. First, for relativistic 
neutrinos, there is no mixture if all mixing matrices are the same: 

rr[(e"'™""'')2] = 1, 

if Ur = Vr = Vl = Ul- We sec also that the size of the mixture depends on the shape of the chosen 
mixing matrices. For some mixing matrices, the mixture can be large and of the order of 0(e^, Vr^Vl)- 
Both sources of state mixtures, helicity mixing and mass mixing, can have comparable effects. For 
example, the dotted-dashed line in Fig. 4, where wc take Ur = Ul Vl = Vr, indicates that the 
departure of Tr[(£i"''""'*'')^] from unity is maximal. Here, we parameterize the V matrices in the same 
way as the standard Ul matrix and take sin^6i2 = 0.87, sm^^23 = 0.09, sm^^is = 0.23 and Scp = 0.5. 
If, instead, we take Vl = Vr = Ul 7^ Ur (the continuous line in Fig. 4), the mixture caused by mixing 
between different states of neutrino masses is marginal. The neutrino mixing in the scalar sector has a 
decisive role. 



3 Conclusions 

We have constructed the density matrix for neutrinos created in any production process. Numerical 
calculations have been performed for the elementary reaction e~ + u ^ i>i + d. The density matrix 
has been calculated in the total Center of Mass system. Beyond this, we have shown that the Lorentz 
transformation to Laboratory system (the rest frame of a detector) has negligible effects on the density 
matrix for experimentally interesting neutrinos. For neutrinos with energies close to their masses, their 
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states must be described by a mixed density matrix. In this region, the effect depends on the hghtest 
neutrino mass, the mass hierarchy, and the particles momenta uncertainties given in wave packets. It was 
shown that, for rclativistic neutrinos, mixing of the states depends on the mechanism of the production 
process. For the vSM, where only a Left-Handed vector current is in effect, neutrinos are produced in 
a pure state with extreme accuracy. In this way, the assumption about so-called flavour states as in 
the Maki-Nakagawa-Sakata-Pontccorvo mixture of the mass states, is justified. If NP interactions also 
influence the neutrino production process, the assumption about purity of neutrino production states is 
not legitimate. Departure from pure states is caused by two things: the possible production of neutrinos 
in two hclicity states caused, for example, by Right-Handcd currents or by scalar interactions, and 
by their mixing given by different mixing matrices in the various sectors of the production mechanism. 
Separate density matrices for the helicity subsystem and mass subsystem were constructed. Both of these 
effects have been investigated numerically and analj^tically in the case of relativistic and non-relativistic 
neutrinos. 
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